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Abstract 



Correct identification of the Bravais lattice of a crystal is an important early step in struc- 
ture solution. Niggli reduction is a commonly used technique. We investigate the boundary 
l ■ polytopes of the Niggli-reduced cone in the six-dimensional space G* by organized random 

probing of regions near 1-, 2-, 3-, 4-, 5-, 6-, 7- and 8-fold boundary polytope intersections. 
We limit our consideration of valid boundary polytopes to those avoiding the mathematically 
r^ ' interesting but crystallographically impossible cases of zero length cell edges. Combinations 

Oh| of boundary polytopes without a valid intersection or with an intersection that would force 

a cell edge to zero or without neighboring probe points are eliminated. 216 boundary poly- 
topes are found. There are 15 5-D boundary polytopes of the full G* Niggli cone, 53 4-D 
Cy ^ boundary polytopes resulting from intersections of pairs of the 15 5-D boundary polytopes, 

79 3-D boundary polytopes resulting from 2-fold, 3-fold and 4-fold intersections of the 15 
5-D boundary polytopes, 55 2-D boundary polytopes resulting from 2-fold, 3-fold, 4-fold and 
higher intersections of the 15 5-D boundary polytopes, 14 1-D boundary polytopes resulting 
from 3-fold and higher intersections of the 15 5-D boundary polytopes. The classification of 
the boundary polytopes into 5-, 4-, 3-, 2- and 1-dimensional boundary polytopes corresponds 
well to the Niggli classification and suggests other possible symmetries. 

All of the primitive lattice types can be represented as combinations of the 15 5-D bound- 
|/-\ \ ary polytopes. All of the non-primitive lattice types can be represented as combinations of 

^' . the 15 5-D boundary polytopes and of the 7 special-position subspaces of the 5-D boundary 

f^i^ ' polytopes. This study provides a new, simpler and arguably more intuitive basis set for 

the classification of lattice characters and helps to illuminate some of the complexities in 
Bravais lattice identification. The classification is intended to help in organizing database 
searches and in understanding which lattice symmetries are "close" to a given experimentally 
determined cell. 



c^ ! 1 Introduction 



The wide-spread use of Niggli reduction in crystallography implies that it should be thoroughly 
understood. Robust identification of Bravais lattices and lookup of unit cell parameters in 
databases would be improved if a successful embedding of the Niggli space could be achieved. 
An essential preliminary step towards doing such an embedding is the mapping of the boundaries 
of the space of Niggli-reduced cells. 

Failure to correctly identify the Bravais lattice of a crystal can compromise subsequent least- 
squares calculations or even the solution of a structure. There are two commonly used ways to 
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obtain a unique representative of the infinite number of cells that may be used to generate a 
given lattice: Niggli reduction [TS] and Delaunay reduction [7]. We follow the conventions of 
the International Tables for Crystallography [B] in basing this paper on Niggli reduction, and we 
will recast the discussion in terms of Delaunay reduction in a future study. Given a precisely 
determined reduced cell, the lattice symmetry may be unambiguously inferred. In addition, 
reduced cells are useful in searching for sub- and super-cells, indexing of powder patterns and in 
twinned crystal studies. Without a suitable metric and a clear understanding of the geometry of 
the boundaries of the space, searches of lattices in the neighborhood of a given cell must either be 
excessively broad and produce many false positives, or, if made tighter, risk missing important 
candidates, especially in the vicinity of 90 degree angles. 

As noted in [3] the concept of a reduced cell is strongly related to the concept of a reduced 
ternary quadratic form |17| [18] . Reduced cells have become an important computational tool in 
crystallography (just as reduced quadratic forms are an important tool in computational number 
theory) but much of the literature focuses on an essentially qualitative approach, taking a reduced 
cell as an absolute indicator of symmetry and not considering the impact of experimental error. 
Consideration of experimental error requires use of a distance metric dealing with reduced cells. 
Such metrics have been used in a 6-dimensional approach |20| . another 6-dimensional approach 
[2], a B-matrix approach [T3] and a "distortion index" approach [Tl]. In addition, because the 
processes of both Niggli and Delaunay reduction can produce large discontinuities in reduced 
cell parameters from small changes in the lattice, an effective use of a metric must allow for 
such discontinuities, either by a combinatorial search or by a metric-preserving embedding in a 
higher dimensional space that removes the discontinuities. Until now, crystallographic software 
appears either to have used a demanding combinatorial approach or simply to have given up on 
a doing complete search. The purpose of this paper is to take the necessary first steps towards 
adding an embedding to the existing combinatorial approaches by clearly mapping the boundary 
discontinuities and their related transformations. 

There have been multiple investigations of such boundaries, albeit without a consideration of 
experimental error in most cases. See, for example, Gruber [5] [lOj for a review and an approach 
using a 5-dimensional space based on the metric tensor. Gruber's 1997 approach partitions the 
space of reduced cells into 127 disjoint components (genera), based on 67 1- to 4-dimensional 
"hyperfaces" further subdivided into 227 hyperfaces in order to achieve a common partitioning 
applicable to both Niggli and Delaunay reduction. Unfortunately, Gruber's reduction to 5 di- 
mensions, and any purely topological approach without a metric that allows error propagation 
from the experimental data, makes it difficult or impossible to do a full perturbation analysis of 
the impact of experimental errors. 

Therefore in this investigation we return to the full 6-diniensional space [2] , G , of unit cells 
based on the metric tensor and use algebraic and Monte Carlo techniques to explore the "natural" 
5-, 4-, 3-, 2- and 1-dimensional boundary polytopes of the 6-dimensional polytope of Niggli- 
reduced cells. The two techniques are mutual supportive. The Monte Carlo technique is especially 
important in clearly delineating the 5-dimensional boundaries, which can be seen to align with 
the Niggli reduction steps. The lower dimensional boundaries are an algebraic consequence of the 
5-dimensional boundaries. The lower dimensional boundaries derived algebraically are confirmed 
by the Monte Carlo technique which helps to identify unpopulated boundaries and boundaries 
that drop to lower dimension due to glancing intersections with multiple other boundaries. 

There are 15 5-dimensional boundary polytopes and a total of 216 boundary polytopes. In 
this approach, all the boundary polytopes are on the surface of the 6-dimensional Niggli-reduced 
polytope. Identification of these boundary polytopes and the reduction transformations involved 
in crossing them is an essential step in either a combinatorial error analysis or in embedding 
G into a higher dimensional space for an analytic error analysis, and it is useful for database 
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searches. 

The 15 5-D boundary polytopes give the complete shape of the space of Niggh-reduced ceUs. 
All of the primitive lattice types can be represented as combinations of the 15 5-D boundary 
polytopes. All of the non-primitive lattice types can be represented as combinations of the 15 
5-D boundary polytopes and of the 7 special-position subspaces of the 5-D boundary polytopes. 
By confining our attention to just the Niggli reduction, the result is a simpler classification than 
Gruber's with more direct applicability to an embedding and database searching. The embedding 
and database searches are the subjects of following studies. 



2 The space G^ 

G^ is a reformulation of the crystallographic metric tensor and the "Niggli matrix" (itself a 
reformulation of the metric tensor) [2]. A vector g in G® is defined as: 

t ^ ( "^•"^, t-t, ~t-~^, 2t-~^, 2~^-~^, 2~^-~t\ 

= { \\ltr, ||tf, lll^f, 2\\t\\\\l^\\cos{a), 
2||^||||^||cos(/3), 2\\t\\\\t\\cos{j)) 

= {91, 92, 53, ff4, 55, 56) =.9{1,2,3},5{4,5,6} 

where a , b , c arc the unit cell edge vectors, and "•" indicates the dot product. The unit cell 
is chosen to be primitive. 

The notation 5/4.5^6} is used to denote the elements (54,55,56) from the full G vector. 

3 The NiggU Conditions 

The Niggli- reduced cell of a lattice is a unique choice from among the infinite number of alternate 
cells that generate the same lattice [TS]. A Buerger-reduced cell for a given lattice is any cell 
that generates that lattice, chosen such that no other cell has shorter cell edges [S]. Even after 
allowing for the equivalence of cells in which the directions of axes are reversed or axes of the 
same length are exchanged, there can be up to five alternate Buerger-reduced cells for the same 
lattice [HI ■ The Niggli conditions allow the selection of a unique reduced cell for a given lattice 
from among the alternate Buerger reduced cells for that lattice. 

Niggli reduction consists of converting the original cell to a primitive one and then alternately 
applying two operations: conversion to standard presentation and reduction |15j [H Table 1]. 
The convention for meeting the combined Buerger and Niggli conditions is based on increasingly 
restrictive layers of constraints: 

If 51 < 52 < 53, I54I < 52, I55I < 51, |56| < 51 and either 5(4,5, 6} > or 5(4,5,6} < then we 
have a Niggli-reduced cell, and we are done. 

The remaining conditions are imposed when any of the above inequalities becomes an equality 
or the elements of 5(4,5.6} sue not consistently all strictly positive or are not consistently all less 
than or equal to zero. 

The full set of combined Buerger and Niggli conditions in addition to those for the cell edge 
lengths being minimal is: 

require < 51 s^ 52 ^ 53 
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if 51 = 52, then require [54! =% I55I (3.1) 

if 52 = 53, then require jgsl ^ \ge\ (3.2) 

require {34 > and ^5 > and 176 > 0} 

or require {(74 ^ and 55 ^ and g6 ^ 0} 

require I54I ^52 

require \g5\ ^ gi 

require {gel ^ 5i 

require 53 ^ 51 + 52 + 53 + 54 + 55 + 56 

if 54 = 52, then require ge < 2^5 (3.3) 

if 55 = 51 , then require gg ^ 234 
if 56 = 51, then require 55 < 2^4 
if 54 — —52 , then require 56 = 
if 55 = —51 , then require 56 = 

if 56 = -51, then require 55 = (3.4) 

if 53 = 51 + 52 + 53 + 54 +55+56, then require 2gi + 2^5 + gg ^ 

The G® transformations associated with each of these steps are enumerated in [5] . Apphcation 
of these operations must be repeated until all are satisfied. 

4 Notation and Boundary Polytopes 

We define the manifold of Niggli reduced cells in G^ as N. 

The interior of N, int{N), is defined as the set of ri £ N such that there exists r > such 
that for ah ~f € G^ such that W^ -lt\\ <r,^ G N. 

The closure of iV, cl{N), is defined as the set of 5 G G^ such that for all r > there exists 
it e iVsuch that ||^ - 7t|| ^ r. 

The boundary of AT, d{]Sl), is defined as the set of points in d(7V) not in int{N). 

The boundary of N is created by the linear constraints of Niggli reduction and therefore can 
be decomposed into the union of "polytopes," i.e. flat facets with straight edges. 

We distinguish the primary boundary polytopes from their edges, which are also polytopes, 
by "dimension" which is the number of vectors in a basis for the interior of the polytope. 

AT is a six-dimensional (6-D) polytope. iV is a double-ended cone-like region going through 
the origin to infinity in both directions. The boundary polytopes are flat facets created by 
the intersections of hyperplanes through the origin. The boundary polytopes are, of course, of 
lower-dimension than N. Therefore any randomly selected vector in G has a vanishingly small 
probability of occupying any particular 5-dimensional ("5-D") boundary polytope, and it has an 
even lower probability of occupying one of the lower-dimensional boundary polytopes resulting 
from the intersections of 5-D boundary polytopes. In addition, some Niggli-reduced boundary 
polytopes are not, themselves, in N. We call those boundary polytopes "open", i.e. while 
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there are Niggli-reduced cells near that boundary, some or all of the points on those boundary 
polytopes are not themselves Niggli-reduced. 

Our task is to identify the 5-D boundary polytopes that give JVits shape. Those 5-D bound- 
aries and the transforms involved in crossing them generate all the rest of the structure. However, 
in order to understand the shape of a given 5-D boundary polytope, we need the 4-D edges that 
bound it. In order to understand the shape of a given 4-D boundary polytope, we need the 3-D 
edges that bound it. In order to understand the shape of a given 3-D boundary polytope, we need 
the 2-D edges that bound it. In order to understand the shape of a given 2-D boundary polytope, 
we need the 1-D edges that bound it. From this classification we gain a better understanding 
of the relationships among Bravais lattice types, and, perhaps more importantly, this provides 
essential information needed to organize computations. Hosoya [llj addressed a different, but 
related, classification. We discuss the relationship between boundary polytope identification, lat- 
tice types and Hosoya's approach in section 1121 Hosoya recognized the complexity of boundary 
identifications for iVand introduced the use of Monte Carlo searching to clarify the relationships. 
We use a Monte Carlo search in section [5j 

The reduction steps convert a non-reduced cell into one that has at least one edge shorter than 
the starting edges, and other steps in the case of equality convert a non-reduced cell into a cell 
that is more orthogonal than the starting cell. These operations are accomplished by choosing 
a face- or body-diagonal to replace one of the cell edges. The conditions added to remove the 
ambiguities in the case of equalities allow for a unique choice of Niggli cell in all cases but thereby 
create complex boundary conditions [2]. For example, the cell edge equalities in equations (13. ip 
and (|3.2p create boundary polytopes across which elements of both gsi^2,3} Etnd of (7{4,5,6} £^re 
exchanged, while equations (|3.3|) through (|3.4p create boundary polytopes across which edges 
are exchanged for face-diagonals. 

A boundary polytope does not necessarily consist entirely of Niggli-reduced cells, but must 
be "near" Niggli-reduced cells and must be "near" non-Niggli-reduced cells as well. We define a 
boundary polytope, F C N as a subset of N for which there is an associated matrix Mr such 
that for all points ^ £ F there exists S > such that for all nn ^ iV where \\nn — 7 || < 
S, Mr ( 7 ) G N. Each boundary polytope is a portion of the boundary for which there is a 
single transformation matrix that maps the nearby non-Niggli-reduced cells to Niggli-reduced 
cells. This is not necessarily a mapping back to the starting point, not even to a point near the 
starting point, not even for points on the boundary polytope itself. We look for point-by-point 
invariance in defining special-position subspaces in section [S] 

5 The Monte Carlo Search for Boundary Polytopes 

In theory, one could follow an entirely algebraic process to identify all boundary polytopes. In 
practice, using a computer to explore N and provide information on the boundary polytopes is 
helpful. 

In order to reduce the computational burden and deal with the open boundary polytopes, 
rather than randomly generating points in G , we generate random short lines in G and look 
for cases in which the boundary of N must have been crossed because one end of the line is 
Niggli-reduced, while the other end is not. The process for the initial search for 5-D boundary 
polytopes is given in Table [T] 

The search for the boundary polytopes resulting from the process in Table [1] produces trans- 
formations in the course of Niggli reduction (step 6). We sort the Niggli-reduced G^ vectors by 
the associated transformation matrix from step 6. A high population for a given matrix indi- 
cates a significant volume of G with access to the associated boundary polytope, implying a 
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Table 1: Initial process to locate 5-D boundary polytopes 

1. Generate a random vector in G . 

2. Confirm that the vector represents a proper unit cell (for instance, the sum of the 
interaxial angles must be less than 360 degrees). 

3. Niggli-reduce the vector. 

4. Randomly perturb the vector resulting from step 3. 

5. Confirm that the perturbed vector represents a proper unit cell. 

6. Niggli-reduce the vector from step 5, accumulating the total transformation matrix 
from the vector of step 3 to the new reduced vector. 

7. If the transformation is the unit matrix, then the perturbed vector was not near 
the boundary of N. Discard the trial. Otherwise, proceed. 

8. If the transformation matrix has been discovered before, increment a counter for 
its occurrences. Otherwise, add it to the list of discovered transformations. 

9. Repeat many times, starting at step 1. 



5-D boundary polytope (See Fig. [T]and Fig. ^. A lower population for a given matrix possibly 
implies a lower-dimensional boundary polytope, a common "edge" resulting from the intersec- 
tion of multiple 5-D boundary polytopes. With enough probes, many of the lower dimensional 
boundary polytopes can be discovered, but only with difficulty because the information on those 
lower-dimensional boundary polytopes is swamped in a sea of data about the 5-D boundary 
polytopes. 

Efficient discovery of all these intersections of multiple 5-D boundary polytopes requires the 
use of "projectors" . A projector for a subspace Q in G^ is a symmetric matrix, P, that when 
applied to any point g in G is such that Pg is the closest point in Q to g . A projector has 
the nice property that PP = P, i.e. that it acts like the identity matrix on the subspace Q. 

The projector for a given boundary polytope can be discovered by examining the set of 
Niggli-reduced G vectors associated with the transform matrix for that boundary polytope. The 
examination can either be a simple inspection of the list of vectors or can be done algorithmically 
by use of a singular value decomposition (SVD) calculation [U [T^l HH] . It should be noted that 
the projector projects onto the hyperplane associated with the polytope and may project a point 
from G^, or even from N itself, outside of N. 

Let us first derive projectors by simple inspection. A search for points near a boundary 
produces a list of vectors that can be examined for the conditions that are met. For instance, 
the trial vectors: 

-2.73956 -1.78806 
-0.26549 -4.79911 
1.57819 1.85157 



4.41605 


53.21164 


53.3171 


-9.85206 


4.95245 


106.2402 


106.5968 


-72.3608 


5.62821 


98.26772 


98.36612 


24.37056 



are seen to meet the condition (72 — .93- 
The trial vectors: 



4.85822 


9.79018 


40.14963 


-3.6758 


-0.01092 


-2.18456 


4.89205 


21.22063 


75.92303 


6.01777 


0.05752 


2.78514 


5.03365 


26.32789 


46.84058 


-17.2646 


-0.01252 


-2.03757 



are seen to meet the condition g^ — 0. 
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Figure 1: Illustration of Monte Carlo search for boundary polytopes in G^. The hatched line 
represents a boundary polytope, e.g. g2 =54- The unhatched side of the boundary (iV) consists 
of Niggli-reduced cells. In this case the hatched side is divided into two regions, NNl {e.g. 
95 ^ ffe) and NN2 {e.g. g^ < go) for which different reduction matrices are required to get a 
Niggli-reduced ceU, e.g. Mq for NNl and My for NN2. ''NN' stands for not Niggli-reduced. VI 
is a randomly generated probe point in region N for which a random short line reaches across the 
boundary to reach region NNl, so the starting point is associated with Mg. V2, VA and V5 are 
randomly generated probe points for which the random short line remains in the Niggli-reduced 
region N. Therefore V2, VA and V5 are discarded. V3 is a randomly generated probe point in 
region N for which a random short line reaches across the boundary to reach region NN2, so 
the starting point is associated with A/y. Because of this difference in reduction matrices, the 
boundary polytope is treated as consisting of two distinct boundary polytopes, in this case, cases 
6 and 7. 






^- V4 
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Figure 2: Counts of points found near various boundary polytopes in 100 million trials, organized 
in declining order of counts, showing the most populated 23 of the 92 boundary polytopes found 
in this run. This is a run with no filtering for any particular boundary with the counts shown 
on a logarithmic scale. Note the precipitous drop of nearly 2 orders of magnitude after the 
first 15 boundary polytopes. This drop confirms that those 15 boundary polytopes are the 5-D 
boundary polytopes and that there is a vanishingly small probability of there being any other 
5-D boundary polytopes of the Niggli-reduced cells. 
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In each case the necessary projector is obvious (see projectors P2 and P3, below). In other 
cases, such as the body-diagonal boundary polytope, or many of the lower-dimensional boundary 
polytopes involving multiple face-diagonals, simple inspection is challenging. Given a sufhcient 
number of vectors, the projector may be recovered algorithmically from such a list of vectors, 
rather than by inspection. If a singular value decomposition is done on a list of vectors near to 
and covering a single 5-D boundary polytope, the vector corresponding to the smallest singular 
value is the unit normal to that boundary polytope, and the projector onto that polytope is 
simply the identity minus the projector onto the line of the normal. The projector onto the 
line of the unit normal is generated by forming a matrix. A, whose first row is the unit normal 
and whose remaining rows are set to zero. Then A"^A is the projector onto the line of that unit 
normal. 

In the general case of a lower-dimensional boundary polytope, the projector onto the boundary 
polytope is the identity minus the projector Q onto the hyperplane f2 orthogonal to the boundary 
polytope. In that case a singular value decomposition on the list of vectors near to and covering 
the lower-dimensional boundary polytope will have small singular values for the vectors spanning 
fi. The projector Q onto f2 is generated by forming a matrix A with initial rows consisting of 
the vectors corresponding to the small singular values and the remaining rows set to zero. Then 
A^A is the projector onto f2 and / — A'^A is the projector onto the boundary polytope. 

Once we have found the projectors for the 5-D boundary polytopes, the projectors for their 
intersections may be generated from their products. The product of two different projectors 
is not, itself, likely to be a projector, but repeated squaring of that product rapidly converges 
to the correct projector [T*. Given these projectors, the process for the remaining boundary 
polytopes is given in Table [21 This process is similar to the initial process, but three new steps 
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Table 2: Process to locate lower-dimensional boundary polytopes. 

0. Read in the projector for the desired boundary to be probed. 

1. Generate a random vector in G . 

2. Confirm that the vector represents a proper unit cell (for instance, the sum of the 
interaxial angles must be less than 360 degrees). 

2A. Project the vector confirmed in step 2 onto the boundary. 
2B. Confirm that the projected vector is still a valid unit celL 

3. Niggli-reduce the vector. 

4. Randomly perturb the vector resulting from step 3. 

5. Confirm that the perturbed vector represents a proper unit cell. 

6. Niggli-reduce the vector from step 5, accumulating the total transformation matrix 
from the vector of step 3 to the new reduced vector. 

7. If the transformation is the unit matrix, then the perturbed vector was not near 
the boundary of N. Discard the trial. Otherwise, proceed. 

8. If the transformation matrix has been discovered before, increment a counter for 
its occurrences. Otherwise, add it to the list of discovered transformations. 

9. Repeat many times, starting at step 1. 



have been added: In step the projector for a particular boundary polytope is read in. In step 
2A the vector confirmed by step 2 is projected onto the hyperplane containing the boundary 
polytope, and step 2B verifies that the projected vector is still a valid unit cell. When doing a 
random perturbation of a vector already projected onto the hyperplane containing a boundary 
polytope while looking for intersections with that boundary polytope, it is sufficient to use a 
random vector projected onto that the hyperplane containing the boundary polytope. However, 
that risks missing some very low dimension cases. That search produced 215 distinct boundary 
polytopes. The search was then redone, stepping back into the Niggli cone after each projection, 
and using a full 6-D spherical random search to ensure catching any nearby boundary polytopes. 
That search added one more distinct boundary polytope, for a total of 216 boundary polytopes. 

The dimension of a boundary polytope can be determined by computing the number of 
eigenvalues equal to 1 of the projector onto the hyperplane containing the boundary polytope, 
and having distinct boundary projectors is a sufficient condition for two boundary polytopes 
to be distinct. However, distinct projectors is a necessary, but not sufficient, condition for two 
boundary polytopes to be distinct, because crossing one bounding hyperplane in two different 
places can require two different transformation matrices to reduce the result. The reduction 
transformation matrices themselves are needed to disambiguate cases with the identical bounding 
hyperplane projector. 

In the course of the Monte Carlo investigation of the boundary polytopes resulting from 
combinations of 5-D boundary polytopes, we track the relative populations as an indicator of 
a consistent assignment of dimension. The candidates for assignment to a particular dimen- 
sion are sorted by population, and the list is cut off on a precipitous drop of more than an 
order of magnitude. For the remaining populations, the mean population /i and estimated stan- 
dard deviation a are computed. For each population r, a Z-score, (r — /i)/(T, is computed (see 



http://en.wikipedia.org/wiki/Standard-SCorep. Z-scores of less than -1 suggest the possible need 



for further study of the boundary polytopes involved. A combination of the population analysis 
resulting from the Monte Carlo search, algebraic analysis of constraints and computation of the 
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eigenvalues of projectors allows the identification of the dimensions of the boundary polytopes, 
helps to identify distinct combinations of intersections of 5-D boundary polytopes that represent 
the same lower dimensional boundary polytopes and confirms the critical initial identification of 
the 15 5-D boundary polytopes identified purely by the Monte Carlo search. 

We limit our consideration of valid boundary polytopes to those avoiding the mathematically 
interesting but crystallographically impossible cases of zero length cell edges. Combinations of 
boundary polytopes without a valid intersection or with an intersection that would force any 
of ff{i,2,3} to zero or that did not have neighboring Niggli-reduced probe points are eliminated. 
574 combinations of 1 through 8 intersecting 5-D boundary polytopes were not degenerate. 
Many combinations represent equivalent boundary polytopes. There are 216 distinct boundary 
polytopes. There are 15 5-D boundary polytopes of the full G^ Niggli cone, 53 4-D boundary 
polytopes resulting from intersections of pairs of the 15 5-D boundary polytopes, 79 3-D boundary 
polytopes resulting from the 2-fold and higher intersections of the 15 5-D boundary polytopes, 
55 2-D boundary polytopes resulting from 2-fold and higher intersections of the 15 5-D boundary 
polytopes, 14 1-D boundary polytopes resulting from 3-fold and higher intersections of the 15 
5-D boundary polytopes. The ability of the intersection of only 2 5-D boundary polytopes to 
produce 4-D, 3-D or 2-D boundary polytopes results from the additional constraints imposed by 
Niggli reduction. 

For the 5-D boundary polytopes, three special cases arise in which the Niggli reduction 
conditions divide a single polytope into two sub-poly topes, the "flat boundary intersections" . In 
each those three cases, instead of two boundary polytopes meeting at an angle, two boundary 
polytopes meet edge on. These three cases are 54 = (?2, <?5 = <7i, and ge = gi- In each case 
the division in the polytope is based on the equality of the other two members of 5(4,5.6} • For 
example, when 94 = 52, the division is along 55 = ge- Therefore, when multiplying projectors, 
instead of multiplying the projector for 54 = 52 by itself when the two half-polytopes are required, 
the second projector is replaced by the projector for the division, i.e. in the example given instead 
of multiplying the projector for 54 = 52 by itself, it is multiplied by the projector for 95 = gg. 

6 Special-Position Subspaces 

In an analysis of symmetry, special positions play an important role. A special position is a 
point invariant under a symmetry transformation, an eigenvector, of eigenvalue 1, of the trans- 
formation. We define a special-position subspace of a boundary polytope, F, as the intersection 
of the eigenspace of eigenvalue 1 of Mr with the boundary polytope. Formally, for a boundary 
polytope, r, with associated matrix, Mr, the special-position subspace A{r, Mr) is defined as 
the set of points 7 S P such that MrCf) = 7 . In the case of boundary polytopes associated 

with a transformation matrix that goes from the all acute -|- -|- -|- case to the all obtuse case 

or vice versa, there cannot be any Niggli-reduced special-position subspace, because the axial 
planes of the 3(4,5,6} subspace are excluded from the all acute + + -I- case. As we will see, while 
the special-position subspaces of the boundary polytopes are not needed in order to classify the 
primitive lattice types, they come into play in classifying the non-primitive lattice types. 

If we have two boundary polytopes Pi , ^"2 C N, we denote the intersection of the closure 
of Fi and the closure of F2 as F1F2. Note that intersection is a commutative operation, i.e. 
FiF2 = F2Fi. 

Because we have restricted the boundary polytope in this paper to have only one associated 
matrix, we use the notation F for A{F, Mr). In general, an infinite number of higher dimensional 
polytopes will intersect F in F. We distinguish such a higher dimensional polytope as F' . Thus 
FF' = F. 
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6 SPECIAL-POSITION SUBSPACES 



Table 3: Fifteen 5-D boundary polytopes on Niggli- reduced cells in G . The special-position 
subspaces are identified by the conditions to be added to the conditions that define the bound- 
ary polytope itself. For a given boundary polytope F, the column "Condition" gives the G^ 
constraints (prior to closure) of the boundary polytope. When taken with the "Special-Position 
Subspace" constraint in the last column, the result is the entirety of the special-position sub- 
space r. The "Special- Position Subspace" constraint by itself is F' . Boundary polytopes 1 and 

2 apply in both the all acute (-1- -I- -I-) and all obtuse ( ) branches of the Niggli-reduced 

cone. Boundary polytopes 8, B, E and F are restricted to the all obtuse ( ) branch of the 

Niggli-reduced cone, N. Boundary polytopes 6, 7, 9, A, C and D are restricted to the all acute 
(-|- -I- -I-) branch of TV. While the boundary polytopes 3, 4 and 5 are boundaries of both the all 

acute (+ + +) and all obtuse ( ) branches, the common special position subspace of those 

polytopes is just 54 = 55 = 56 = which is part of the ( ) branch. 



Class 


Case 


Condition 


Special- Position 


Subspace 


Equal cell edges 


1 


51 =52 


.94 = 55 






2 


52 = 53 


.95 = 56 




Ninety degrees 


3 


.94 = 


55 = 56 = 






4 


55 = 


54 = 56 = 






5 


.96 = 


54 = 55 = 




Face diagonal 


6 


52 = 54 and 55 > ge 


(none) 






7 


52 = 54 and 55 < ge 


55 = 56/2 






8 


.92 = -.94 


(none) 






9 


51 = 55 and 54 > ge 


(none) 






A 


51 = 55 and 54 < ge 


54 = 56/2 






B 


51 = -55 


(none) 






C 


51 = 56 and 54 > 55 


(none) 






D 


51 = 56 and 54 < 55 


.94 = .95/2 






E 


51 = -56 


(none) 




Body diagonal 


F 


51 + 52 + 53 + 54 + 55 + 56 = 53 


51 - .92 - .94 + .95 = 


= 
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7 The fifteen 5-D boundary polytopes 

The fifteen 5-dimensional boundary polytopes and their special-position subspaces may be or- 
ganized as shown in Tabled in which we use the hexadecimal digits 1 through F as identifiers. 
For each 5-dimensional boundary polytope F in Table [3] having a non-trivial special-position 
subspace, we designate the particular choice of higher dimensional polytope intersecting F in F 
as F'. See section [7.0.31 for a concrete example. 



7.0.1 Equal cell-edge cases: 

Cases 1 and 2 arise when two cell edges have equal length. The conditions of Niggli reduction 
impose a secondary condition on the associated angles for those two edges that resolves the 
ambiguity in ordering them. For example, incase 1, || a || = || 6 || (gi = .g2) and the Niggli- reduced 
G'* vector U = (gi,g2,5'3,g4,g5,g6) produces the same lattice as V = (g2, gi,53,g5,g4,56), 
which is not Niggli-reduced if g^ and 35 have different values. In going from Niggli-reduced cells 
near U to Niggli-reduced cells near V (e.g. by decreasing g2 slightly), we are crossing a boundary 
polytope with a discontinuity in each of (74 and 175 . We may represent the transformation that 
takes U into V at the first boundary polytope by the matrix Mi that maps U into V and the 
projector Pi that maps any G vector into the gi = (72 boundary polytope. 



Case 1: 51 =52, 



~ct -i' b , 



^^ 



-4 



Ml = 



/O10000\ 

100000 

10 

000010 

000100 

V 1 y 

Gr subspace: (r, r, s, i, u, v) 



Pi = 



/ .5 


.5 











\ 


.5 


.5 




















1 




















1 




















1 





V 














1/ 



Similarly, for case 2, || & || = || c || (32 = .93), 55 and g^ are exchanged, yielding 
Case 2: 52 = ff3, b ^ 1^ , ~^ 



t 



Mo 



/ 1 \ 

001000 
010000 

10 

000001 
yooooioy 

G subspace: (r, s, s, i, u, v) 



P-?. 



( 1 





Vo 




.5 
.5 







\ 

5 

5 

10 

10 

1 y 



The remaining cell-edge case in which || a | 

\t\\ = \\-t 



II c II ( (71 = (73), is only considered under Niggli 
reduction when ||7?|| = || b || and || 6 || = ||~cr|| which is a combination of case 1 and case 2. This 
requires two simultaneous 5-dimcnsional constraints, thereby making g\ = 33 a 4-dimensional 
rather than a 5-dimensional case. 

The special-position subspaces 1 and 2 are obtained by adding the constraints 1' : {54 = 175} 
and 2': {175 — gg}, respectively. 
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7.0.2 Ninety-degree cases: 

Cases 3, 4 and 5 arise when a reduced-cell angle is ninety degrees. In those cases, the remaining 
cell angles both can be replaced by their supplements. This changes the sign of <?{4,5,6}- 



34 ==0, it ^ 


-7? 










/ 1 





\ 




/ 1 


\ 


10 










1 





10 
1 










P3 = 


1 










-1 










1 


\ 





-W 




i^ 


1/ 



M, = 



G subspace: (r, s,t,0, —u, —v) 



Case 4: 


55 = 0, b 


^-6 












/ 1 





\ 




/ 1 


\ 




1 










10 





Mi = 


1 
0- 



-1 






Pi = 


10 
1 











1 















i^ 





-w 




i^ 


l) 



G subspace: (r, s, t, — u, 0, —v) 



Case 5 



9e 


= 


T 


^- 


- c 




/ 1 














\ 





1 




















1 




















-1 




















-1 





Vo 














1/ 



M5 = 



G subspace: (r, s, t, 



P^ = 



/ 1 \ 

10 

10 

10 

10 

\ooooooy 



-v,0) 



In each ninety-degree case, the special-position subspace consists of 3, 4, 5 : {(74 — g^ ~ gg ~ 
0}, i.e. the primitive orthorhombic case, and we take 3' : {55 = gg = 0}, 4' : {54 — g^ = 0}, 
5': {54 =.95 = 0}. 

7.0.3 Face-diagonal cases: 

Cases 6 through E arc all face-diagonal cases, in which a cell edge is equal in length to a face- 
diagonal. Some complexity arises in the analysis because, unlike Delaunay reduction, Niggli 
reduction permits non-obtuse angles. We can always change the sign of any two elements of 
9{4 5 6) by changing the direction of the cell edge involved with those two elements. For example, 
if we transform a to — a then, then while gi remains unaffected, the signs of each of 35 == 2 a-c 
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.^t 



and qq — 2 a ■ b will change. Thus we can transform a cell having three acute angles to a cell 
having one acute angle and two obtuse angles, and we can transform a cell having three obtuse 
angles to a cell having one obtuse angle and two acute angles. The complete list of possible sign 
changes in 5/4.5^6} by changing the directions of axes are: 



{+ + +} 



{--+} 



{- + -} 



{+--} 



{- 



-} 



{+ + -} 



{+-+} 



{- + +} 



Unless one of the angles is ninety degrees (which introduces a zero into .g{4,5.6}), we cannot 
ordinarily transform a Niggli-reduced cell with all-acute angles to one with all-obtuse angles by 
changing the directions of axes, nor can we transform a Niggli-reduced cell with all-obtuse angles 
to one with all-acute angles by changing the directions of the axes. Note that changing the 
direction of all three axes has no effect because all the sign changes cancel. 

The face-diagonal cases do include cases in which transformations from, for example, -I- -I- -I- 
to do occur. Let us look in detail at cases 6 and 7, .92 ~ 54- 



tt 



52 

52 -54 
52 - 54 + 53 



-> 



-> 



-^ -^ 



lt-7^lp 



ll^ 



^1 



54 



53 

ll^ll 



-> 



to b 



-^ 



Thus transforming c' 

gs are, of course, unchanged and 



will not change the cell edge lengths. In this case, gi, 52, and 



54 

9', 
9', 



1^) 



2't-Ct 
2-^-Ct --^ 



■i^ 



-> 



2d- b 



56 



2t^ 



-^ 



-> 



2 b- c = 252 - 54 = 54 



= 21^- b - 21^-1^ ^ gQ- g^ 



This shows that a single element of 5{4,5,6}, 55, will change sign depending on the sign of 56 — 55- 
Cases 6 and 7 cannot start from the all-obtuse case because g^ = 32 and because 32 niust be 
nonnegative. Starting from an all acute case, + + +, we will remain in the all acute case if g^ 
is greater than or equal to 55 but go to one having one obtuse angle (not reduced) if g^ is less 

than (75. We then change to having all obtuse angles, , by reversing the direction of b. The 

resulting matrices in these face-diagonal cases are: 



Case 6 



52 = 


= 54, 55 


^56, 


b 


-^-; 


7, ^- 


^ b 


-^ 




/ 1 











\ 




( 1 





\ 





1 
















.5 


.5 






1 1 

-2 


-1 
1 










P6 = 






1 
.5 



.5 











-1 


1 










1 


yo 











-W 




^0 





1/ 



Mfi 



G subspace: {r,s,t^s,u-\- v,v) 
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Case 7: 52 = 94, 9b < 56, "^ 



-> 



1^ 



M7 = 



/ 1 





\ 




/ 1 





\ 


1 













.5 


.5 


1 1 
2 


-1 
-1 






P7 = 






1 
.5 



.5 








-1 1 










1 


i^ 





l) 




^0 





1/ 



G subspace: (r, s, t, s,M, u + w) 



Mr 



52 = -54, a , 


^- 


-7?, 


Y ^ b 


+ d 






/ 1 





\ 




( 1 





\ 


10 













.5 


-.5 


111 
2 1 










Ps = 






1 
-.5 



.5 





-1 


-1 










1 


\ 





-w 




lo 





1 7 



G subspacc: (r, s, t, — s, 



Mo = 



= 55, 9i 


^56, 


1^- 


-^-"fe", ^- 


->^- 


-1^ 






1 








\ 




/ .5 








.5 \ 


1 
















1 








1 1 





-1 


-1 






1 


^9 = 






1 





1 






2 





1 







.5 








.5 











-w 




I 








l) 



V 



G subspacc: {r,s,t,u-\- v,r^u) 



Case A: 51 =55, 94 < 96, 



Ma 



/ 1 





\ 




/ .5 


.5 \ 


1 










10 





1 1 





-1 


-1 
1 


Pa = 


10 
1 






2 





-1 




.5000 


.5 


\ 





1) 




i^ 


1) 



G subspacc: (r, s,t,u,r,u + v) 



Case B: gi 



-55, 



t 



-1^, 



1^ 



7? 



^ 
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Mb = 



/ 1 

1 

1 1 





1 



2 

V 

G subspace: (r, s, t, 



1 








0-10-1 







Pb = 



0-1/ 



.5 











.5 


\ 





1 




















1 




















1 








.5 











.5 




















1/ 



-¥ 



/ 1 














1 


1 

















1 

















-1 


1 














-1 


1-2 















Case C: gi =96,94 > 95, b 



Mc: = 



G subspace: (r, s,i,M + w, w, r) 



Case D: gi = ge, 94 < gb, b 



Mr 



G subspace: {r,s,t,u,u + v,r) 



-1 




1/ 



-^ 



, b -^ 


b --^ 






/ .5 


.5 \ 




10 





Pc = 


10 
10 








1 







i^ .5 


W 



^ 



/ 1 














1 


1 

















1 

















-1 


1 














1 


^2 















\ 


/ 


1 








Pd = 







w 


\ 



5 

10 

10 

10 

1 



.5\ 








.5 / 



Case E: 


91 = -96, 


1;^^ 


1^ + t 


, ^^ 


-^ 






/ 1 





\ 




/ .50000 


-.5 \ 




1 1 





1 




10 





Me = 


1 





-1 - 



-1 


Pe = 


10 
10 











- 


-1 




1 







i^ 2 





1) 




i^ -.5 


.5/ 



G subspace: (r, s, t, — u, — w, — r) 

The special-position subspaces of the face-diagonal boundary polytopes 6, 8, 9, B, C and E 
are empty because such a special position would require a common point in the all acute + + + 

and all obtuse cases that only meet at the axial planes of the (7{4.5,6} subspace, which are 

excluded from the all acute -|--l--|- case. For cases 7, A and D there are non-trivial special-position 
subspaces. An invariant point in case 7 would have to satisfy g^ = gg — 55 or g^ ~ .96/2- Thus 
we define 7' : {(75 = 56/2} and similarly define A' : {g^ = .96/2} and D' : {g^ = .95/2}. 
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8 THE 4-D BOUNDARY POLYTOPES 



Table 4: The 53 4-D boundary polytopes. The boundary conditions use the symbol notation 
of Table |3] for the 5-D boundary polytopes. The bounding 4-D polytopes of each of the 15 5-D 
boundary polytopes can be read off from this table by going down the appropriate column to 
the main diagonal and then across. For example, the edges of the 3 polytope are 13, 23, 34, 35, 
3A, 3B, 3D and 3E. 



92 = 93 


.94 = 


55 = 


56=0 


a-face-diagonals 


b- face-diagonals 


c-face-diagonals 


body 


12 


13 
23 


14 
24 


15 
25 


16 17 18 
26 27 28 


19 lA IB 
29 2A 2B 


IC ID IE 
2C 2D 2E 


IF 
2F 






34 


35 

45 


47 48 
56 58 


3A 3B 

59 5B 


3D 3E 

4C 4E 








67 


69 


7C 


8F 








9A 


9C 

AD 

BE 


BF 










CD 


EF 



7.0.4 Body-diagonal case: 



^11 



There is only one 5-dimensional body-diagonal case, || a + 6 + c || = || c || : 



CaseF: 51+52+53 + 54 + 55+56=53, ~^ ^ -{ 



( 



Mf = 



V 



1 


1 

-2 






1 
1 1 

-2 







1 
-1 







1 
-1 



-1 -1 



Pf = 



G subspace: (r, s, i, — u, — u, 



1 / 

S -f M -I- w) 



It 



( -8 
-.2 



-.2 

-.2 

V--2 



^) 



-•2\ 
-.2 


-.2 
-.2 

.8/ 



In order to have a special-position subspace in case F, in addition to 51+52+53+54+.95+56 = 
^3, we need (from the fourth and fifth rows of Mp), g^ = —2g2 — 54 — 56 and 175 — —2gi — 55 — 9e 
from which we have 2g2 + 2g4 = — gg = 2.gi + 2g<^ from which we take F' : {gi —52 — 54 + 55 = 0}. 
This is equivalent to || a + c || = || 6 + c ||, i.e. that the shorter b-face-diagonal is the same 
length as the shorter a-face-diagonal. 



8 The 4-D boundary polytopes 

The 4-D boundary polytopes are created by the intersection of 2 5-D boundary polytopes. Certain 
intersections are degenerate. For example, cases 8, B, E and F are restricted to the branch 
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of the boundary of the NiggU cone while cases 6, 7, 9, A, C and D are restricted to the + + + 
branch. More subtly, cases 6 and 7 require (72 =34, maximizing 54, forcing 54 ^ g^ and 54 ^ ge 
which would conflict with cases A and D without cases 9 and C, respectively. Similarly, cases 
9 and A maximize (75, which would conflict with case 7 without case 6, and cases C and D 
maximize ge which would conflict with case 6 without case 7. Thus cases 6A, 7A, 6D, 7D, 79, 
7A, 6C, 6D actually are the lower dimension cases 69A, 79A, 6CD, 7CD, 79A, 79A, 6CD and 
6CD, respectively. This process can result in 3-D or even 2-D boundary polytopes from the 
intersection of two 5-D boundary polytopes (see below) . After excluding the cases that involve 
any of f;{i^2.3} = 0, there are fifty-five 4-D cases as shown in TableS) The relative populations for 
all the 2D boundary polytopes except 26, 28, 2 A and 2D have Z-scores above -1. The Z-scores 
for those 4 cases range from -1.1 down to -1.9. 

The edges of the 5-D polytopes can be read directly from Table H) For example the 6 polytope 
is bounded by 16, 26, 56, 67 and 69 and the F polytope is bounded by IF, 2F, 8F, BF and EF. 
It is important to note that the polytopes 1, 2, 3, 4 and 5 extend into the boundaries of both 

the + + + and branches of N. Even though the polytopes 3, 4 and 5 do not contain any 

valid Niggii reduced + + + cells, they are part of both branches of d{N) even for + -f- -I-. 

9 The 3-D boundary polytopes 

The 3-D boundary polytopes are created by the intersection of 3 5-D boundary polytopes. In 
some cases the boundary polytope is better represented by a 4-fold intersection. The bound- 
ary polytope 34CD is equivalent to 34C and 34D, 359A is equivalent to 359 and 35A, 4567 is 
equivalent to 456 and 457, 679C is equivalent to 69C and 79C, and 9ACD is equivalent to 9 AD 
and ACD. These are "flat boundary intersection" cases in which one side of the flat boundary 
intersection implies the other. On the other hand 126 and 127, 12A and 129, 12C and 12D, 2AD 
and 29C, 69C and 79C are distinct rather than equivalent pairs of flat boundary intersections. 
Six 2-fold intersections of 5-D boundary polytopes (6A, 6C, 79, 7D, 9D and AC) result in 3-D 
boundary polytopes, rather than in 4-D boundary polytopes. In each case both boundary poly- 
topes have mismatched partners from "flat boundary intersections" . Let us examine the 6A case 
in detail, elaborating on the discussion of the 4-D boundary polytopes above. 

Cases 6 and A are 52 = 54 and g^ ^ ge and gi — 55 and 5(4 < ge, respectively. For 
intersections, we use the closures of the boundary polytopes, so we have the closure of A as 
51 = 55 and (74 ^ 36- The Niggii cone itself imposes the additional restrictions g^ < gi and 
5i ^52, but from 6, 52 = 34 and from the closure of A, 54 ^ ge^ so we have 

52 = 34 ^ 36 ^ 31 => 32 ^ 31 
and from the Niggii reduction conditions 

31 =^32 

thus gi = (72 and 34 = ge, meaning that, in addition to satisfying the constraints of case A, we 
also satisfy the constraints of case 9, 31 = 35 and 34 ^ ge- Thus case 6 A is actually case 69 A, 
producing a true 3-D boundary polytope from the intersection of 2 5-D boundary polytopes 
due to the additional constraints of Niggii reduction. As we will see in the discussion of the 
2-D boundary polytopes, the Niggii reduction constraints can result in shedding one or more 
degrees of freedom, allowing some 2-fold intersections of 5-D boundary polytopes to result in 
2-D boundary polytopes. 

After excluding the cases that involve any of 3{i,2,3} ~ 0, there are seventy-nine cases as 
shown in Table \5\ The relative populations for all the 3-D boundary polytopes have Z-scores 
greater than -1.1. 
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Table 5: The 79 3-D boundary polytopes. In some cases the most natural presentation of a given 
3-D polytope is a 4-fold intersection. In each of those case an equivalent 3-fold intersection is 
given in parentheses immediately below the 4-fold. 



54 = 55 = .96 = 


a- face-diagonals 


b- face-diagonals 


c-face-diagonals 


body 


123 124 125 
134 135 

145 


126 127 128 

147 148 
156 158 
167 


129 12A 12B 
13A 13B 

159 15B 
169 

19A 


12C 12D 12E 

13D 13E 

14C 14E 

17C 

IAD 

ICD 


12F 

18F 

IBF 
IFF 


234 235 
245 


247 248 
256 258 
267 


23A 23B 

259 25B 
269 

29A 


23D 23E 
24C 24E 

27C 

29C 

2AD 

2BE 
2CD 


28F 

2BF 
2EF 


345 




359A 35B 

(359) 


34CD 34E 

(34C) 

3AD 

3BE 






4567 
(456) 

458 




47C 48EF 
(48E) 








569 

58BF 

(58B) 










679C 
(679) 












9ACD 

(9AC) 












BFF 
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For completeness, if one wished to include the boundary polytope with 51 = it would be 
considered in the 3-D polytopes. In d{N), gi = forces Qb = 96 — 0, ^2 — 0, leaving only 3 
degrees of freedom at d{N). Note that §2 = is of even lower dimension, 1, because §2 — 
forces 54 = as weU as gi = g^ — g^ — 0, leaving only 1 degree of freedom (g^). ^3 = is just 
the origin. 

10 The 2-D boundary polytopes 

The 2-D boundary polytopes are, in general, created by the intersection of 4 5-D boundary poly- 
topes, but several well-populated 4-fold intersections result in 3-D boundary polytopes rather 
than 2-D boundary polytopes. Several 4-fold intersections are most naturally presented as higher 
multiplicity intersections, and in some cases the intersection of as few as 2 5-D boundary poly- 
topes is sufficient to create a 2-D boundary polytope. After excluding the cases that involve any 
of 3(1. 2, 3} — 0, there are fifty-five cases as shown in Table |6l The combination of 7 boundary 
polytopes 1679ACD is the hexagonal rhombohedral hR, lattice character 9, Roof/Niggli sym- 
bol 49B, subspace {r,r,s,r,r,r). Alternatively, lattice character 9 can be viewed as any of 81 
other intersections, including 2 2-folds (6D, 7A), 18 3-folds (179, 16A, 16C, 16D, 17A, 17D, 19D, 
lAC, 67A, 67D, 69D, 6AC, 6AD, 6CD, 79A, 79D, 7AC, 7AD), 33 4-folds (1679 167A, 167C, 
167D, 169A, 169C, 169D, 16AC, 16AD, 16CD, 179A, 179C, 179D, 17AC, 17AD, 17CD, 19AC, 
19AD, 19CD, lACD, 679A, 679D, 67AC, 67AD, 67CD, 69AC, 69AD, 69CD, 6ACD, 79AC, 79AD, 
79CD, 7ACD), 21 5-folds (1679A, 1679C, 1679D, 167AC, 167AD, 167CD, 169AC, 169AD, 169CD, 
16ACD, 179AC, 179AD, 179CD, 17ACD, 19ACD, 679AC, 679AD, 679CD, 67ACD, 69ACD, 
79ACD) and 7 6-folds (1679AC, 1679AD, 1679CD, 167ACD, 169ACD, 179ACD, 679ACD) and 
is a very highly populated 2-D boundary polytope. If we exclude this case, the remaining fifty- 
four cases have Z-scores ranging from -1.36 (for 1456) to 0.84 (for 123E) to 2.75 (for 29ACD). 

Let us consider how one of the 2-folds, 6D, results in only 2-degrees of freedom. Cases 6 and 
D are (72 — ff4 and 55 ^ ge and gi = ge and 54 < 55 , respectively. The closure of D is (71 = ge 
and (74 ^ (?5, and Niggli reduction requires ffe ^ 5i ^ .92, from which it follows that 

56 = 51 5% 52 = 54 < 35 < 51 
and therefore 

56 = 51 = 52 = .g4 = .g5 
creating the subspace (r, r, s,r, r, r), i.e. 2 degrees of freedom. 

11 The 1-D boundary polytopes 

There are 14 distinct 1-D boundary polytopes, with many equivalent presentations. The most 
complex situation is best presented as an 8-fold intersection, 12679ACD, i.e. <?i = 52 = 53 = 
9i = 9b = 56j which is the face-centered cubic {r,r,r,r,r,r). There are 81 other equivalent 
presentations of the face-centered cubic, inherited from the 7-fold intersection hexagonal rhom- 
bohedral discussed above by adding case 2 to each of those presentations, thereby providing 2 
3-folds for the face centered cubic (26D and 27A). The remaining 13 1-D boundary polytopes 
are 12345 (r,r,r, 0, 0, 0), 1234CD (r,r,r, 0, 0,r) (equivalent to 1234C, 1234D, 123CD, 124CD), 
1234E (r, r, r, 0, 0, -r), 12359A (r, r, r, 0, r, 0) (equivalent to 12359, 1235A, 1239A, 1259A), 1235B 
(r,r,r,0,-r,0), 123AD (r, r, r, 0,r,r), 123BEF (r, r, r, 0, -r, -r) (equivalent to 123BE, 123BF, 
123EF, 12BEF, 23BEF), 124567 (r,r,r,r, 0, 0) (equivalent to 12456, 12457, 12467, 12567), 12458 
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1 1 THE 1-D BOUNDARY POLYTOPES 



Table 6: The 55 2-D boundary polytopes 



35 = 56 = 


a-face-diagonals 


b-face-diagonals 


c-facc-diagonals 


body 


1234 1235 
1245 


1247 1248 
1256 1258 
1267 


123A 123B 

1259 125B 
1269 

129A 


123D 123E 
124C 124E 

127C 

12AD 
12CD 


128F 

12BF 
12EF 


1345 




134CD 

(134C) 

1359A 135B 
(1359) 


134E 

13 AD 13BEF 
(13BE) 






14567 1458 
(1456) 




147C 

148EF 
(148E) 








1569 

158BF 
(158B) 












1679ACD 
(6D, 7A) 




2345 




2359A 235B 
(2359) 


234CD 234E 

(234C) 

23AD 

23BE 






24567 2458 
(2456) 




247C 

248EF 
(248E) 








2569 

258BF 

(258B) 










2679C 
(2679) 










29ACD 

(29AC) 














2BEF 
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12 RELATIONSHIP BETWEEN BOUNDARY POLYTOPES AND LATTICE-TYPE 



(r,r,r,-r, 0,0), 1247C (r, r, r, r, 0,r), 1248EF (r, r, r, -r, 0, -r) (equivalent to 1248E, 1248F, 
124EF, 128EF), 12569 (r, r, r, r, r, 0), 1258BF (r,r,r, -r, -r,0) (equivalent to 1258B, 1258F, 
125BF, 128BF). 

These 14 1-D boundary polytopes of JV correspond exactly to the 14 vertices of the hyperpoly- 
hedra in [HI Table 1] for which none of (7{i,2,3} are zero. The 14 that match are a confirmation 
of the completeness of this analysis. Due to the distortion introduced by projection, the re- 
jection of the cases for which any of 5{i,2,3} are zero is important in preserving the metric for 
incommensurate edges near the origin. 

If we exclude the highly populated face-centered cubic, the Z-scores for the relative popula- 
tions of the remaining 13 1-D boundary polytopes range from more than -1.1 for 12569 to 1.44 
for 123BEF. 

12 Relationship between boundary polytopes and lattice- 
type 

"Lattice characters" provide a finer-grained division of lattice type than the 14 Bravais lattice 
types [SI International Tables, Volume A]. In order to understand the relationship among between 
the 216 G^ boundary polytopes and the 44 lattice characters in the International Tables, we use 
combinations of the 15 5-D boundary polytopes and of the special-position subspaces of those 
polytopes. There are multiple alternate representations of some of the lattice characters. We 
discuss some of them below. 

We refer to Roof's redrawn Niggli figure identifiers [TB] as "Roof/Niggli symbols." We 
may associate the Roof/Niggli symbols, lattice characters and Bravais lattice types with the 
indicated subspaces of G^ and combinations of boundary polytopes and special conditions as 
shown in Tables [7] and [8] The triclinic lattice characters 31 and 44 are not included because no 
boundary polytopes are needed for the triclinic case as they fill the Niggli-reduced cone. 

The primitive Bravais lattice types have a simple relationship to the boundary polytopes. 
The primitive cubic, which has one degree of freedom as a G subspace, is the intersection of 
five 5-D boundary polytopes. The primitive tetragonal and primitive hexagonal lattice types 
each have two degrees of freedom as G^ subspaces and each is the intersection of four 5-D 
boundary polytopes. The primitive orthorhombic lattice type has three degrees of freedom as a 
G subspace and is the intersection of three 5-D boundary polytopes. The primitive monoclinic 
lattice types each have four degrees of freedom as G^ subspaces, and each is the intersection of 
two 5-D boundary polytopes. 

The combination of 8 boundary polytopes 12679ACD (equivalent to the 3-folds 26D and 27A) 
is the face-centered cubic cF, lattice character 1, Roof/Niggli symbol 44C, subspace (r, r, r, r, r, r) 
[2]. Alternatively, cF can be viewed as any of 127, 12 A or 12D and of several other intersections. 
As one would expect from the large number of intersecting boundary polytopes, this is a very 
complex region of G and will be the subject of a later paper. 

The 1-D combinations of 6 boundary polytopes 1234CD (equivalent to the 5-folds 1234C, 
1234D, 123CD, 124CD), 12359A (equivalent to the 5-folds 12359, 1235A, 1239A, 1259A), 123BEF 
(equivalent to the 5-folds 123BE, 123BF, 123EF, 12BEF, 23BEF), 124567 (equivalent to the 5- 
folds 12456, 12457, 12467, 12567), 1248EF (equivalent to the 5-folds 1248E, 1248F, 124EF, 
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1 2 RELATIONSHIP BETWEEN BO UNDAKY POLYTOPES AND LATTICE- TYPE 



Table 7: Roof/Niggli symbol, International Tables (IT) lattice character, Bravais lattice type, 



G subspace [2 , G boundary polytope. 



Roof/ 

Niggli 

Symbol 


IT 

Lattice 
Char 


Bravais 

Lattice 

Type 


G' 
Subspace 


G' 
Boundary 
Polytope 


44A 


3 


cP 


(r,r,r, 0,0,0) 


12345 = 123=124=125 


44C 


1 


cF 


(r, r, r, r, r, r) 


12679ACD 


44B 


5 


cl 


(r,r,r,-2r/3,-2r/3,-2r/3) 


12F2'F' = 12F 


45A 
45B 


11 
21 


tP 
tP 


(r,r,s, 0,0,0) 
(r,s, 5,0,0,0) 


1345 = 13 = 14 = 15 

2345 = 23 = 24 = 25 


45D 
45D 

45C 
45E 


6 
7 

15 
18 


tl 
tl 
tl 
tl 


(r, r, r, — r + s, — r + s, —2s) 

(r, r, r, -2s, -r + s, -r + s) 

(r, r, s, — r, — r, 0) 

(r, s, s,r/2,r,r) 


12FF' = 12F 

12F2' = 12F 

158BF 

2ADA' = 2AD 


48A 
48B 


12 
22 


hP 
hP 


(r,r,s, 0,0, — r) 
(r,s,s,-s,0,0) 


134E 
2458 


49C 
49D 
49B 

49E 


2 
4 
9 
24 


hR 
hR 
hR 
hR 


(r, r, r, s, s, s) 

{r,r,r, — s, — s, — s) 

(r, r, s, r, r, r) 

(r, s, s,-s + r/3,-2r/3,-2r/3) 


121'2' = 12 

121'2' = 12 

1679ACD 

2i^2'F' = 2F 


50C 


32 


oP 


(r,s,t, 0,0,0) 


345 = 3 = 4 = 5 


50D 
50E 
50A 
50B 
50F 


13 
23 
36 

38 
40 


oC 
oC 
oC 
oC 
oC 


(r,r,s,0,0, -t) 
(r,s,s,-t, 0,0) 
(r,s,i,0, -r,0) 
(r,s,i,0,0,-r) 
(r,s,i,-s, 0,0) 


134 
245 
35B 
34E 

458 


51A 
51B 


16 
26 


oF 
oF 


(r,r,s,-i, -i, -2r + 2t) 
(r, s,i,r/2,r, r) 


IFl' = IF 
ADA' = AB 


52A 
52B 
52C 


8 

19 
42 


ol 
ol 
ol 


(r, r, r, — s, — t, — 2r + s + t) 

(»■, s, s, i, r, r) 

{r,s,t, —s, —r,0) 


12F 

29C = 2AD 

58BF 
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1 2 RELATIONSHIP BETWEEN BO UNDAKY POLYTOPES AND LATTICE- TYPE 



Table 8: Roof/Niggli symbol, International Tables (IT) lattice character, Bravais lattice type, 
G^ subspace, G^ boundary polytope, continued 



Roof/ 


IT 


Bravais 


G' 


G" 


Niggli 


Lattice 


Lattice 


Subspace 


Boundary 


Symbol 


Char 


Type 




Polytope 


53A 


33 


mP 


(r,s,t,0,-M,0) 


35 


53B 


35 


mP 


{r,s,t,-u,0,0) 


45 


53C 


34 


mP 


(r,s,t,0,0, -u) 


34 


57B 


17 


ml 


(r, r, s, —t, ~u, — 2r + t + u) 


IF 


57C 


27 


ml 


{r,s,t,u,r,r) 


9C = AD 


57A 


43 


ml 


(r, s, t, —s + u,—r + u, —2m) 


FF' = F 


55A 


10 


mC 


(r, r, s, t, t, u) 


11' = 1 


55A 


14 


mC 


{r,r,s,~t,-t,-u) 


11' = 1 


55B 


20 


mC 


{r,s,s,t,u,u) 


22' = 2 


55B 


25 


mC 


(r, s, s, —t, —u, —u) 


22' = 2 


56A 


28 


mC 


{r,s,t,u,r,2u) 


AA' = i 


56C 


29 


mC 


{r,s,t,u,2u,r) 


DD'=D 


56B 


30 


mC 


{r,s,t,s,u,2u) 


77' = 7 


54C 


37 


mC 


{r,s,t, — M, — r, 0) 


5B 


54A 


39 


mC 


{r,s,t, — M,0, — r) 


4E 


54B 


41 


mC 


{r,s,t, — s, — u, 0) 


58 



24 



13 SUMMARY AND CONCLUSIONS 



Table 9: The Hosoya boundary polytopes of the three cubic lattice types and the primitive 



hexagonal in terms of G 



Hosoyai 


2.91 = 55 + 36 


Hosoya2 


2^2 == 54 + 56 


Hosoya^ 


253 = 54 + 55 


Hosoya4 


54 — 55 , the special-position subspace of case 1 


Hosoya^ 


55 = 56j the special-position subspace of case 2 


HosoyaQ 


54 ~ 55/2, the special-position subspace of case D 



128EF), 1258BF (equivalent to the 5-folds 1258B, 1258F, 125BF, 128BF) form the subspaces 

{r,r,r,0,0,r) 
(r, r, r, 0, r, 0) 
(r, r, r, 0, — r, —r) 
(r, r, r, r, 0, 0) 
(r, r, r, — r, 0, — r) 
(r, r, r, — r, — r, 0) 

of which none are Niggli- reduced and arc therefore a set of open boundary polytopes. 

Most of the 216 G^ boundary polytopes are non-Niggli-reduced open boundary polytopes 
of the Niggli region. Therefore only two of the 5-fold boundary polytopes, five of the 4-fold 
boundary polytopes, eight of the 3-fold boundary polytopes and eleven of the 2-fold boundary 
polytopes correspond directly to lattice characters. None of the single 5-D boundary polytopes 
corresponds to lattice characters. Both the non-Niggli-reduced open boundary polytopes and 
the 5-D boundary polytopes are suggestive of symmetries to be investigated in the future. 

We are working from the boundary polytopes, looking for the resulting symmetries. Hosoya 
[Tl] started, instead, from the three highest symmetry lattice types (the three cubics) and added 
the lower symmetry primitive hexagonal lattice type, and inferred boundary polytopes from 
the symmetries, having to treat open boundary polytopes as if they were Niggli- reduced. The 
resulting six boundary polytopes of the three cubics and the primitive hexagonal restated in 
terms of G^ are given in Table [HI 

13 Summary and Conclusions 

The wide-spread use of Niggli reduction in crystallography implies that it should be thoroughly 
understood. Robust identification of Bravais lattices and lookup of unit cell parameters in 
databases would be improved if a successful embedding of the Niggli space could be achieved. 
We have investigated and enumerated the several kinds of boundary polytopes on the Niggli 
cone and enumerated the transformations and projectors specific to each. While other, related, 
work has often involved 3-D sections, our work natively addresses the boundary polytopes and 
transformations in the Euclidean space G . The single point of view and the simple linear algebra 
involved makes the presentation more consistent. 

Some unexpected complexities have been encountered, such as the occurrence of one boundary 
polytope that is the intersection of eight 5-D boundary polytopes, a fruitful area for further 
investigation. 
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